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Abstract 

We study the effect of wave function orthogonality in the rela- 
tivistic treatment of the nucleon removal reactions (7,p) and {e,e'p). 
The continuum wave function describing the outgoing nucleon is made 
orthogonal to the relevant bound states using the Gram-Schmidt pro- 
cedure. This procedure has the advantage of preserving the asymp- 
totic character of the continuum wave function and hence the elastic 
observables are unaffected. The orthogonality effects are found to be 
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negligible for (e, e'p) reactions for missing momenta up to 700 MeV/c. 
This holds true for both parallel and perpendicular kinematics. By 
contrast the orthogonalization of the wave functions appears to have a 
more pronounced effect in the case of (7,]3) reactions. We find that the 
orthogonality effect can be significant in this case particularly for large 
angles. Polarization of the outgoing protons and photon asymmetry 
show more sensitivity than the cross sections. If the orthogonality con- 
dition is imposed solely on this one hole state the effects are usually 
smaller. 
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1 Introduction 



Nucleoli removal reactions play a central role in clarifying important aspects 
of the nature of nuclear states, particularly those related to the shell struc- 
ture of nuclei. There is currently renewed interest in these reactions when 
the probe is electromagnetic as more sophisticated electron beam machines 
are now available. The new generation of possible experiments opens the 
door to more precise and elaborate data that will help us understand the 
intricacies of the reaction mechanism and look at possible changes in the 
structure of nucleons and other hadrons as they move within the nuclear 
environment. Spin dependent observables are expected to play a key role in 
these investigations. 

Most earlier analyses of reactions such as the quasielastic (e, e'p) and 
(7,p) reactions have been carried out using the non-relativistic distorted 
wave impulse approximation (NRDWIA) . In this approach the nucleon 
wave functions needed for the calculation of the reaction amplitude are so- 
lutions of the Schrodinger equation, and the current operators are obtained 
from the free probe-nucleon relativistic amplitude via some non-relativistic 
reduction scheme 0. Recently however it has become clear that a fully rel- 
ativistic approach to the description of these reactions is more appropriate 
0, 1^, |, ^ H, [1^. This is evidenced by the improvement in the values 
of the spectroscopic factors obtained using this approach and indications of 



improvements to recoil polarization [0, |T^ as well as cross sections at large 
missing momenta |]l3i . 

Because of the expectation that more precise and expanded data will 
become available shortly, there has been concentrated activity aimed at dis- 
cerning details of the theoretical calculations of the amplitudes. The interplay 
between relativistic and non-relativistic approaches has been studied by sev- 
eral authors, and recently aspects of gauge dependence as well as coupling 



effects have been investigated [14, 15 



The purpose of the present paper is to address one detail of the calcula- 
tions of the amplitudes which is usually overlooked: namely, both relativistic 
and non-relativistic DWIA versions of the amplitude for a variety of reac- 
tions, including (e, e'p) and (7,p) (which are the focus of the present study) 
suffer from an orthogonality defect. These distorted wave models utilize 
wave functions for bound and continuum states, which are solutions of equa- 
tions containing different hamiltonians; the bound states are generated from 
real potentials, while the continuum states are generated from complex en- 
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ergy dependent optical potentials. The character of the latter potentials is 
dictated by the need to account for nucleon elastic scattering observables. 
The wave functions are therefore not orthogonal, in contradiction to what a 
fundamental theory would require. 

The impact of the nonorthogonality defect has been discussed in the past 
and ways to remedy it have been suggested |116l O, |18|, EOl pT], p2l |23 



The majority of these discussions were done in a non-relativistic framework 
but many of the remedies suggested have never been utilized in main stream 
calculations. Noble Jl^ attempted to recast the amplitude for photoabsorp- 



tion into a form that takes into account the orthogonality of the nucleon wave 
functions. Difficulties with this approach have been discussed by Eisenberg, 
Noble and Weber |TH]. Celenza and Shakin [O suggested that a modified 



hamiltonian for the continuum states can be constructed in such a manner 
as to guarantee both orthogonality with the bound states, and a fit to the 
proton elastic scattering data. This proposition was never tested in practice. 



Other possibilities discussed by Boffi et al. pT| , yfJl involved the introduc- 
tion of effective operators, but it turned out that these were too difficult to 
calculate. 

For the two reactions under consideration, however, a more practical ap- 
proach was identified and used to assess the role of non-orthogonality on the 
non-relativistic amplitudes. It was suggested by Boffi et al. ||2T|, ^ as well 
as by Ciofi Degli Atti et al. ||20| that the continuum states can be made 



orthogonal to the bound states by applying a Gram-Schmidt orthogonaliza- 
tion procedure. In the present work we use this same method to examine 
the importance of orthogonalization in the relativistic approach to nucleon 
knockout reactions. 

In the following text we give a brief outline of the relativistic calculations 
for nucleon knockout induced by the quasi elastic scattering of electrons or 
the absorption of real photons. We then introduce the method of orthogonal- 
ization of the nucleon wave functions and calculate the resulting corrections 
to the reaction amplitude. A discussion of the impact of these changes on 
the reaction observables is given in section 3. We conclude in section 4. 

2 The Relativistic Calculations 

The relativistic calculations of the amplitude, in the one photon exchange 
model for the (e, e'p) process, are discussed in Refs. 0, ^, H, |^, |^ . Relativistic 
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calculations for (7,^) reactions are discussed in Refs. |T^. It is sufficient 
for the purpose of the present discussion to focus on the nuclear matrix 
element portion of the amplitude, which can be written in a unified form for 
reactions involving either virtual or real photons. In the notation of Ref. 
the relativistic expression for the nuclear matrix element A^^^^^, leading to 
a specific final state of the residual nucleus for both the (7,p) and (e, e'p) 
reactions can be written as 



^\T^J■MB 



J d^r )] ^ (p, r) Tp^j^^MB (r) exp {tq ■ r) , (1) 



where Mb and fi are the spin projections of the bound and continuum nucle- 
ons, respectively, and p is the 3-momentum of the latter. The 4x4 matrix 
r^, operating on the nucleon spinors, is given in detail for the (e, e'p) re- 
action, in Eq. (2.8) of reference [Q. For the (7,p) reaction, is written 
explicitly in Ref. |2|. 

In relativistic DWIA (RDWIA) calculations of the amplitude, is cho- 
sen as the Dirac continuum wave function in the elastic channel. It is gen- 
erated by using a complex nucleon-nucleus optical potential, and is required 
to give an acceptable account of the elastic data. The spinor \E'j^^Mb, on the 
other hand, is a solution of the bound state Dirac equation, which contains 
real scalar and vector potentials. It is clear then that the nucleon wave func- 
tions used in the RDWIA are not orthogonal since they are solutions of the 
Dirac equation for two different hamiltonians. It is the lack of orthogonality 
of the spinors \I'^ and "^jg^MB that is the subject of the present study. 



2.1 Orthogonal States 

The orthogonal continuum wave function is obtained via the Gram-Schmidt 
method. This involves subtracting a multiple of the bound state wave func- 
tions for each of the occupied levels of the target nucleus. The approach has 
the advantage of being easily incorporated within the RDWIA approach as a 
means of gaining some insight into the possible corrections arising from the 
orthogonality defect. The proposed orthogonal wave function can be written 
in terms of the original distorted wave and the bound state wave functions 
as 

^{+) (p, r) = (p, r) - Y: Pi'-''-'''-^L',,j',,M'B (r), (2) 

L'^^J'bM'b 
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where the sum is over all the occupied bound states (characterized here 
by their angular momentum labels L'^, J'^, M'^). Note that the coefficient, 

^' ^, in the second term of Eq. (0) is required to have the same asymp- 
totic spin projection fi, as the original continuum wave function. 

It was shown by Ciofi Degli Atti et al. [20| that there is a close connection 
between the above form of orthogonalized wave functions and the require- 
ment of antisymmetrization of the final state wave function in reactions ini- 
tiated by the real photons. A crucial advantage of the above prescription is 
that the resulting orthogonal continuum wave functions continue to yield the 
same proper description of elastic scattering provided by the original proton 
continuum wave function. 

The expansion coefficients (3] 
the wave function given in Eq. 
target; this yields 



are determined by requiring that 
is orthogonal to all occupied states of the 



Lb,Jb,Mb l^J 



(3) 



The detailed evaluation of these expansion coefficients requires the stan- 
dard partial wave expansion for the wave functions. For the continuum state 
we write 

1/2 

Aw 



(+) 



2m 



LJM 



Y, 



N 

M-fi 



(L, 1/2; M-/i, /i I J, M) 



3^ 



M 
L 1/2 J 



(4) 



hj (r) 

where the upper and lower component radial functions, fij (r) and Qlj (r), 
are solutions of equations containing complex scalar and vector potentials. 
For the bound state with angular momentum quantum numbers Lb, Jb, Mb 
we can write 

IlbJb i^) 
-WLbJb (^) ^ ■ 



"^JbMb 



-\i Mb 
^Lb 1/2 Jb 



(5) 



Using Eqs. (B) and 



Lb,Jb,Mb 



in Eq. (|^), we get for the expansion coefficients 



2m 



Mb- 



N 



{Lb,1/2;Mb- ^i,I^\Jb,Mb) 



^Lb,Jb- 



(6) 
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where bi^^jg is the overlap integral between the bound and distorted wave 
radial functions: 



ftTf (r) ftj^ (r) + g^Tf (0 £ii (^)l • (7) 



Using Eqs. (|) and (|^) in Eq. (g) allows us to write a partial wave 
expansion for the distorted wave which will be orthogonal to the bound 
states. The result is that the continuum wave function is changed only in 
partial waves with angular momentum quantum numbers corresponding to 
the occupied levels of the target nucleus. The resulting continuum wave 
function is suitable for use in calculations of the nuclear matrix element of 
Eq. (|l]) without further modification, and is straightforward to include in 
existing reaction codes. 



3 Discussion 

To assess the impact of the wave function nonorthogonality in realistic sit- 
uations we have carried out calculations for a number of cases involving 
photonuclear knockout of protons for photon energies above 60 MeV. We 
also extend the study to electron quasielastic scattering on the same nuclei. 
In these calculations we use bound state wave functions which are generated 
using phenomenological Dirac scalar and vector potentials with Woods-Saxon 
shapes. The parameters of the potentials are determined by the requirement 
that the separation energies of the states for which we have experimental 
data are reasonably well reproduced. This binding potential is used to gen- 
erate the entire family of occupied states in the target. The distorted waves 



are generated using the energy dependent potentials of Cooper et al. ||25 
For the ^^O target we use the energy-dependent potentials specific to this 
target (EDAI 016), while for the other nuclei considered we use the energy- 
and mass-dependent set (EDAD Fit 1) These distorted wave functions are 
then made orthogonal to the family of bound states following the procedure 
discussed in section 2. 

The size of the orthogonality effect is linked to the magnitude of the coef- 
ficients b^gjg given in Eq. (0). As is seen from this equation, the coefficients 
are radial overlaps between bound and continuum wave functions. Figure 1 
shows the behavior of the overlap coefficient for the lpi/2 state in ^^O as a 
function of the energy of the continuum nucleon. Part (a) of the figure shows 
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the real (solid) and imaginary (dashed) parts of the coefficient on a stan- 
dard plot, while part (b) shows the absolute value of the real and imaginary 
parts on a semilogarithmic plot. As one might expect these overlap coeffi- 
cients are larger at lower energies and drop quickly as the energy increases. 
They have oscillatory behavior reflecting that particular character of the con- 
tinuum radial wave functions. At low energies (~10 MeV) the magnitude 
of the coefficients can be close to 0.1. Although these are relatively large, 
they are not likely to have large effects in the energy region of interest for 
the reactions under consideration. Note also that the orthogonality require- 
ment affects only a limited number of partial waves (depending on the target 
nucleus) in the continuum wave function: those corresponding to occupied 
states in the target nucleus. 

In the present work, the effects of orthogonality are studied in two cases: 

i) Minimal orthogonality: the continuum wave function is made orthogo- 
nal to only the one-hole state involved in the reaction; 

ii) Full orthogonality: the continuum wave function is made orthogonal to 
all occupied states of the target. 

In the following we give an account of the results obtained for the two reac- 
tions under consideration. 

3.1 {'J,p) Reactions: 

We have investigated the effect of non-orthogonality on the cross section, 
proton polarization and photon asymmetry in (7,p) reactions on four nuclei: 
i^C, ^°Ca and ^°^Pb, for photon energies from 60 to 500 MeV. Two 
representative examples of these calculations are shown in Figs. 2 and 3. 
In Fig. 2 we show the calculations for the reaction ^^O (7,p) at an incident 
photon energy of = 100 MeV. A proton is removed from the lpi/2 level 
and the residual nucleus is left in the ground state. In Fig. 3 we show 
similar calculations for incident photon energy E-y = 196 MeV. To put the 
comparisons in a proper context we also show the experimental cross section 
data for the above reactions |2^, |27[] . 

We note that the effects of minimal orthogonalization are generally smaller 
than those for full orthogonalization. The cross sections (part (a) of the fig- 
ures) are affected only slightly at most angles but the deviations can be large 
at extreme backward angles. The relative increase of the effect at large angles 
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is a reflection of tlie fact, mentioned above, that the orthogonahty require- 
ment modifies only low partial waves. The effects on proton polarization 
(part (b)) and photon asymmetry (part (c)) are noticeable over most of the 
angular range, but are small at the very forward angles. This is an indication 
that these two observables are somewhat more sensitive to the orthogonality 
defect than are the cross sections. 

Figures 4 and 5 shed more light on the role of orthogonality and the 
dependence on the photon energy and the angle of the outgoing protons. 
Here we show the photon energy dependence of the cross section on '^^Ca 
at proton angles of 9p = 90° (Fig. 4) and 9p = 135° (Fig. 5). The cross 
section data shown in both figures are from Leitch et al. [^. We observe 
from Fig. 4 that at 9p = 90° the effects on the direct knockout cross section 
are rather small up to photon energies of 300 MeV (calculations show this 
feature still prevails at higher energies). At 6p = 135° (Fig. 5), we find 
larger effects on the cross section, particularly when full orthogonalization is 
imposed on the wave functions. These results are in qualitative agreement 
with the observations made by Ciofi Degli Atti et al. using the non- 
relativistic approach. Parts (b) and (c) of the figures show the effects on 
proton polarization and photon asymmetry, respectively. Here we note large 
effects on both observables at both angles and over the entire energy range 
(these effects are much smaller at 6p = 45°). It is evident that the spin 
observables are more sensitive to the orthogonality of the wave functions. 

Boffi et al. [^, ^ argued that the effect of non-orthogonality on the 
unfactorized (7,p) cross sections vanishes because of the transverse nature 
of the electromagnetic current. Our present calculations do not support this 
conclusion; as seen in the preceding figures, we do observe some effects when 
using only minimal orthogonalization. The reason for the difference is that 
in Eq. (^), even if one restricts the sum to a single orbital {Lb, Jb) there are 
still contributions from different magnetic substates. As a result some of the 
contributions from orthogonalization do not vanish. This difference can be 
stated as due to the possibility of spin-flip in the outgoing channel (which is 
closely linked to the presence of the spin-orbit interaction) ; such a possibility 
was not taken into account by Boffi et al. 

Earlier estimates of the orthogonality effects in (7,p) reactions were ob- 
tained for non-relativistic calculations by Fink et al. and for relativistic 
calculations by Lotz ^9). The procedure followed was to replace the dis- 
torting potential of the continuum proton with the bound state potential. 
The continuum and bound state hamiltonians are thus identical and hence 
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the wave functions are orthogonal. The dotted curve in Fig. 4 shows the 
calculated observables when this procedure is implemented. The effects are 
seen to be much more exaggerated in this case. It should be noted that one 
important shortcoming of this procedure is that the new distorted wave func- 
tion is no longer able to reproduce the elastic scattering data. This difficulty 
is also present in non-relativistic calculations that use continuum RPA wave 
functions [Q. 



3.2 (e, e'p) Reactions: 

We have also carried out similar calculations for the (e, e'p) reaction. Boffi 
et al. and Ciofi Degli Atti et al. have reported that in the non- 

relativistic calculations of the amplitudes the orthogonalization effects are 
quite small. Our relativistic calculations confirm this conclusion. We have 
carried out calculations for electron energies of 460 MeV and 2.4 GeV and 
found that the spectral functions are affected only slightly by the imposition 
of orthogonality. Our calculations covered missing momenta in the range 
±700 MeV/c. At proton energies of 100 MeV or higher the effects of or- 
thogonality are essentially negligible. At low proton kinetic energies we find 
some orthogonality effects in the proton polarization calculations. Figure 6 
shows calculations for the reaction ^^O (e, e'p) in perpendicular kinematics 
with a proton energy of Tp = 35 MeV. It is evident in this figure that or- 
thogonality effects are small except at the large missing momenta. Figure 7 
shows quasifree electron scattering calculations for the ^^Ca target in parallel 
kinematics, with a proton energy of Tp = 100 MeV. The orthogonality cor- 
rections due to Eq. are negligible, but replacing the distorting potential 
with the binding potential, as was done in Fig. 4, has a large effect. This 
prescription cannot be advocated because the distorted wave is no longer 
able to describe proton elastic scattering. Note that these latter calculations 
cannot describe the electron scattering data very well either. 



4 Conclusion 

In this paper we have investigated the role of wave function orthogonality 
in relativistic calculations of (7,^9) and (e, e'p) reactions. We have used a 



Gram-Schmidt procedure ET], |2^, EO] to orthogonalize the continuum and 



bound-state wave functions of the nucleon. This procedure, besides being 
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closely linked to the antisymmetrization |^D[| of the final state wave function, 
also has the merit of leaving the continuum wave function asymptotically 
intact. This ensures that the wave function continues to provide the same 
good description of the elastic scattering data as the original distorted wave. 
Calculations were carried out for two schemes, one where orthogonality is 
imposed only in relation to the specific hole state involved in the reaction, 
and full orthogonality with all occupied levels. The effects in the latter case 
were found to be larger. 

The present study suggests that the orthogonality effects are not a seri- 
ous problem for (e, e'p) reactions. Calculations over a wide range of missing 
momenta do not show large deviations from the traditional nonorthogonal 
calculations except for the spin dependent observables at large missing mo- 
menta. 

The orthogonality effects are found to be non-negligible in the case of 
(7, p) reactions. Our calculations show that these effects tend to be more pro- 
nounced at large angles of the outgoing protons. The spin observables appear 
to be more sensitive to the orthogonality of the wave functions. This may 
have important ramifications for calculations of the amplitudes for higher 
order processes such as those associated with two-nucleon currents. 
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Figure Captions 



FIG. 1. The overlap coefficients, big^jg, of Eq. (J^) for the lsi/2 state in 
the ^^O target. The binding potentials have Woods-Saxon shapes, while 
the proton optical potentials are the E-dependent optical potential for ^^O 
from Ref. (a) linear plot of real (solid) and imaginary (dashed) parts, 

and, (b) semilogarithmic plot of absolute value of real (solid) and imaginary 
(dashed) parts. 

FIG. 2. Effect of wave function orthogonalization in the ^^O (7, p)^^ A^g.s. re- 
action involving the knockout of a lpi/2 proton. The energy of the incident 
photon is 100 MeV. The binding potentials have Woods-Saxon shapes, while 
the proton optical potentials are the E-dependent optical potentials for ^^O 
from Ref. (a) Cross section, (b) Proton polarization, (c) Photon asym- 
metry. Solid curves — no orthogonality. Long-dashed curves — minimal 
orthogonality. Short-dashed curves — full orthogonality. The data are from 
Ref. p|. 



FIG. 3. Same as Fig. 2, except the energy of the incident photon is 196 
MeV. The data are from Ref. 



g.s 



FIG. 4. Effect of wave function orthogonalization in the '^^Ca{'y,p)'^^ K, 
reaction involving the knockout of a 1^3/2 proton. The angle of the final 
proton is 90°. The binding potentials have Woods-Saxon shapes, while the 
proton optical potentials are E- and A-dependent optical potentials (Fit 1) 
from Ref. |^. (a) Cross section, (b) Proton polarization, (c) Photon asym- 
metry. Solid curves — no orthogonality. Long-dashed curves — minimal 
orthogonality. Short-dashed curves — full orthogonality. Dotted curves — 
continuum wave function calculated using binding potential. The data are 



from Ref. [BS 



FIG. 5. Same as Fig. 4, except the angle of the final proton is 135°. The 
data are from Ref. [BHl. 



FIG. 6. Effect of wave function orthogonalization in the ^^O {e,e'p)^^ N, 



g.s. 



reaction involving the knockout of a lpi/2 proton. The energy of the incident 
electron is 456 MeV, with perpendicular kinematics, and the proton kinetic 
energy is Tp = 35 MeV. The binding potentials have Woods-Saxon shapes, 
while the proton optical potentials are the E-dependent optical potentials 
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for ^^O from Ref. (a) Cross section, (b) Proton polarization. Solid 



curves — no orthogonality. Long-dashed curves — minimal orthogonality. 
Short-dashed curves — full orthogonality. 

FIG. 7. Effect of wave function orthogonalization in the ^°Ca (e, e'p)^^ ^g.s. 
reaction involving the knockout of a ld3/2 proton. The energy of the inci- 
dent electron is 460 MeV, with parallel kinematics, and the proton kinetic 
energy is Tp = 100 MeV. The binding potentials have Woods-Saxon shapes, 
while the proton optical potentials are E- and A-dependent optical poten- 
tials (Fit 1) from Ref. (a) Cross section, (b) Proton polarization. Solid 
curves — no orthogonality. Long-dashed curves — minimal orthogonality. 
Short-dashed curves — full orthogonality. Dotted curves — continuum wave 



function calculated using binding potential. The data are from Ref. [31 . 
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